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-rand  hence  size,  of  low-frequency  arrays  that  our  present  knowledge  of  the  internal  structure  of 
the  deep  ocean  will  permit. 

The  special  objective  of  this  report  is  the  presentation  of  an  algorithm,  the  supporting  theory 
of  which  has  previously  been  published,  which  has  subsequently  been  supported  with  a  series  of 
experimental  investigations  of  the  transverse  horizontal  coherence  properties  of  the  acoustic 
field  following  long-range  propagation  in  the  deep  ocean.  From  this  algorithm  the  gain  of  an  array 
can  be  calculated  provided  it  is  restricted  to  those  instances  where  bottom  and  surface  interaction 
play  a  minor  part  in  the  overall  intensity  of  the  received  field.  <4.  - 
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PREFACE 


This  report  is  a  major  output  of  a  program  in  spatial  properties  of  low-frequency  acoustic  fields  in 
the  deep  ocean.  This  project  was  instituted  at  the  Naval  Research  Laboratory  in  1974  to  provide  a  priori 
estimates  of  the  capabilities  and  limitations  of  large-array  construction  and  performance  due  to  coher¬ 
ence  degradation  from  environmental  causes.  The  project  work  has  emphasized  stochastic  propagation 
measures  of  irregularities  in  the  ocean.  This  approach  has  been  followed  to  provide  probabilistic  predic¬ 
tions  of  the  expected  environmental  limits  to  aperture  designs. 

The  order  in  which  contributing  elements  have  been  dealt  with  in  this  program  has  been  volume 
effects,  bottom  effects,  and  surface  effects.  The  latter  two  categories  are  currently  under  development, 
while  the  first  category  is  dealt  with  in  part  in  this  report  and  is  essentially  closed  at  this  time.  The  or¬ 
dering  of  the  mechanisms  and  their  influence  has  been  guided  by  the  logical  separability  of  volume 
effects  from  those  that  relate  to  the  bottom  and/or  the  surface.  Thus,  the  volume  effects  treated  in  the 
present  work  provide  the  outside  practical  limit  on  resolution,  and  hence  size,  of  low-frequency  arrays 
that  our  present  knowledge  of  the  internal  structure  of  the  deep  ocean  will  permit.  In  other  words,  the 
resolution  limit  caused  by  the  forward  scattering  of  the  acoustic  waves  by  internal  inhomogeneities  in 
the  deep  ocean  can  be  estimated  from  the  results  of  the  present  program.  The  model  provided  in  this 
report  calculates  the  horizontal  component  of  this  resolution  limit. 

The  special  objective  of  this  report  is  the  presentation  of  an  algorithm,  whose  supporting  theory 
has  previously  been  published,  and  which  has  subsequently  been  supported  with  a  series  of  experimen¬ 
tal  investigations  of  the  transverse  horizontal  coherence  properties  of  the  acoustic  field  following  long- 
range  propagation  in  the  deep  ocr  an.  From  this  algorithm  the  gain  of  an  array  can  be  calculated  pro¬ 
vided  it  is  restricted  to  those  instances  where  bottom  and  surface  interaction  play  a  minor  part  in  the 
overall  intensity  of  the  received  field.  The  report  includes  a  detailed  discussion  of  the  mutual  coher¬ 
ence  function  which  characterizes  the  kernel  of  a  probabilistic  array  gain  computation,  followed  by  a 
description  of  the  acoustic  field  in  a  wavenumber  sense.  This  is  elaborated  in  considerable  detail  in  the 
form  of  a  manipulative  program  designed  to  handle  horizontal  line  arrays.  This  development  is 
intended  in  future  works  to  be  enlarged  to  include  the  coherence-degrading  effects  of  vertical  apertures. 
In  its  present  form,  it  provides  a  framework  for  the  inclusion  of  the  forward-scattering  component  ini¬ 
tially  discussed  as  well  as  the  ability  to  include  further  smoothing  factors  that  arise  from  bottom  and 
surface  low-angle  forward  scattering  or  diffusion.  Such  smoothing  factors  are  currently  under  develop¬ 
ment.  Other  program-related  reports  that  extend  the  dimensionality  and  range  of  application  of  these 
calculations  are  listed  as  Refs.  I  through  4. 
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AND  LOW-FREQUENCY-ARRAY  GAIN  LIMITS 
IN  THE  DEEP  OCEAN 


INTRODUCTION 

A  loss  in  the  spatiai  coherence  of  an  acoustic  signal  due  to  random  environmental  variations 
places  an  upper  bound  on  the  performance  of  array  systems.  The  loss  of  coherence  for  totally  refracted 
paths  results  from  a  stochastic  volume  scatter,  the  cause  of  which  can  be  traced  to  internal  temperature 
fluctuations.  For  the  range  of  parameters  of  interest  in  low-frequency  propagation,  these  fluctuations 
can  be  regarded  as  a  direct  result  of  internal  waves.  The  first-order  acoustic-modeling  task  is  to  relate 
the  spatial  coherence  of  the  received  signal  to  the  three-dimensional  spectrum  of  temperature  fluctua¬ 
tions  characteristic  of  internal  waves. 

As  a  consequence  of  the  high  degree  of  anisotropy  of  the  ocean  environment,  which  has  a  hor¬ 
izontal  quasi-lenticular  fine  structure,  loss  in  the  spatial  coherence  of  a  propagating  acoustic  signal 
depends  strongly  on  the  propagation  direction  as  well  as  the  receiver-pair  orientation  relative  to  that  fine 
structure.  Thus,  over  a  vertical  plane  transverse  to  the  direction  of  propagation,  the  rate  of  loss  of 
coherence  for  horizontally  separated  receiving  positions  will  be  much  less  than  that  for  vertically 
separated  receiving  positions.  The  rates  of  loss  along  lines  that  form  an  arbitrary  transverse  angle  with 
the  horizontal  will  fall  between  these  two  extremes.  Further,  the  degree  of  anisotropy  is  such  that  it 
cannot  be  readily  accommodated  by  a  change  of  parameters  in  a  single  acoustic  model.  The  very  struc¬ 
ture  of  the  model  changes  with  the  direction  of  the  line  along  which  the  coherence  is  to  be  estimated  in 
the  transverse  plane.  This,  along  with  the  presence  of  a  sound-speed  field  which  produces  multiple 
paths,  greatly  complicates  the  task  of  estimating  signal  coherence  over  a  vertical  aperture  in  comparison 
to  a  similar  task  for  a  horizontal  aperture.  Reflecting  these  differences  in  the  underlying  physics, 
separate  models  have  been  developed  for  horizontal-transverse  and  for  vertical-transvers ;  spatial- 
coherence  estimation. 

The  model  used  to  estimate  transverse  horizontal  coherence  has  been  designated  COHORT  and  is 
described  extensively  in  the  present  report.  Models  for  the  estimation  of  vertically  transverse  coher¬ 
ence  have  also  been  developed  —  they  are  termed  COVERT  and  CEM.  COVERT  estimates  the 
diffusion  of  a  single  macroray.  and  it  is  most  appropriate  for  acoustic  fields  that  can  be  approximated  by 
a  small  number  of  discrete  paths.  CEM  propagates  the  mutual  coherence  function,  with  scattering,  in  a 
range-depth  plane,  and  thus  it  carries  forward  an  implicit  intensity  field  that  is  smoothed  by  multiple 
low-angle  forward  scattering.  As  might  be  expr  ted,  the  advantage  of  propagation  of  the  coherence 
function  is  purchased  with  additional  computer  a.;  -c-  -  uH  soded  cost. 

The  appropriate  use  of  the  models  described  above  requires  that  the  distinctive  aspects  of  estimat¬ 
ing  coherence  of  an  acoustic  field  be  precisely  delineated.  For  example,  a  conventional  propagation 
model,  incorporating  a  stcpping-parabo'ic-equation  algorithm  (51,  produces  complex  field  points  at  a 
given  range  which  may  be  processed  over  a  depth  interval  to  form  the  complex-conjugate  cross  product 
as  a  function  of  vertical  point  separation.  This  function,  identical  in  form  to  the  mutual-coherence 
function,  contains  only  the  deterministic  part  of  the  acoustic  field;  its  discrete  Fourier  transform  pro¬ 
duces  the  angular  distribution  of  intensity  computed  over  the  selected  depth  interval  of  the  transform. 
This  calculation,  in  othci  words,  models  the  beam  power-output  of  an  array  that  has  receivers  at  the 
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field-point  depths  and  is  equal  in  length  to  the  selected  interval.  The  complex  cross-product  function 
from  such  a  calculation  generally  fluctuates  and  decreases  in  value  with  point  separation,  i.e.,  it  shows 
limited  "coherence"  which  is  entirely  a  result  of  the  interfering  multipath  field.  Actually,  any  set  of  two 
or  more  interfering  waves  produces  "incoherence,”  so  that  an  experimental  measurement  or  a  calcula¬ 
tion  which  accidentally  or  deliberately  includes  multipaths  will  show  "loss  of  coherence."  The  research 
in  the  present  program  has  conceptually  and  practically  separated  this  multipath,  deterministic  incoher¬ 
ence  from  that  generated  by  low-angle  scattering  to  preserve  the  viewpoint  and  potential  value  of  con¬ 
cepts  like  resolution  limit,  which  are  most  often  encountered  in  optical  literature,  or  other  areas  where 
single-path  propagation  is  prominent. 

In  ocean  acoustics,  horizontal  sound-speed  gradients  transverse  to  the  propagation  direction  are 
seldom,  if  ever,  strong  enough  to  produce  horizontal  multipaths.  Thus,  for  a  horizontal  receiving  aper¬ 
ture,  the  single-path  concept  of  resolution  limit  and  transverse  horizontal  coherence  attributable  to  sto¬ 
chastic  mechanisms  may  be  usefully  retained.  The  COHORT  model  calculates  the  coherence  limit  on 
array  signal  gain  for  this  precise  condition.  To  reiterate  for  clarity,  the  aperture  whose  performance 
limit  is  to  be  estimated  is  assumed  to  be  approximately  horizontal  and  transverse  to  the  propagating 
field  whose  horizontal  coherence  is  to  be  calculated.  The  approximate  nature  of  the  orientation  is  just 
that  required  to  constrain  observable  and  deterministic  multipath-induced  coherence  effects  to  some 
required  minimum  level.  The  COVERT  model,  in  a  parallel  conceptual  vein,  estimates  the  single-path 
resolution  limit  vertically  transverse  to  the  propagating  acoustic  field  along  a  macroray.  Thus, 
COHORT  and  COVERT  separately  estimate  the  transverse  horizontal  and  vertical  components  of 
coherence  due  to  forward  scattering,  including  contributions  from  all  the  expected  features  of  an 
ensemble-averaged  picture  of  ocean  acoustic  propagation  dominantly  in  the  volume.  Toe  anisotropy  of 
the  ocean,  with  its  tabular  or  horizontal  lenticular  structure,  plays  a  major  role  in  the  derivation  of  both 
models.  The  principal  inhomogeneous  component  of  the  ocean,  the  vertical  profile,  is  also  included  in 
both  models,  excluding  the  effects  due  to  multipath.  The  Combined  Effects  Model  (CEM)  adds  the  impli¬ 
cit  propagation  of  the  transverse  multipath  structure  to  the  diffusion  induced  by  the  anisotropy  and 
modulated  by  the  inhomogeneity.  It  is  thus  the  most  complex  model  produced  to  date  in  the  present 
program. 

The  COHORT  model,  described  in  detail  in  a  later  section,  is  written  in  a  standard  form  of  FOR¬ 
TRAN  amenable  to  conversion  to  diverse  computer  systems.  It  is  now  in  the  Texas  Instruments 
Advanced  Scientific  Computer  (ASC)  system  at  NRL  and  may  be  accessed  from  remote  terminals.  The 
model  is  also  available  through  NALCON,  a  network  of  interconnected  Navy-laboratory  computers,  for 
remote  operation  or  total  transfer  for  local  operation. 

In  the  remainder  of  this  report,  we  begin  with  a  short  historical  summary  of  propagation  in  ran¬ 
dom  media.  We  continue  with  an  analytic  background  discussion  of  the  coherence  function  that  treats 
the  concepts  described  above  briefly  and  qualitatively.  This  is  followed  by  a  presentation  of  the  model 
that  evaluates  the  horizontal  signal  coherence  and  a  presentation  of  the  wavenumber-space  representa¬ 
tions,  which  are  particularly  useful  for  coherence  estimates.  Beam  patterns  and  array  signal  gain,  the 
desired  final  result,  are  covered  next.  The  next  section  contains  a  detailed  description  of  the  model  and 
its  associated  subroutines.  The  COHORT  model  is  a  driver  program  which  calls  subroutines  in  accor¬ 
dance  with  a  specific  task.  The  program  allows  for  several  optional  entry  and  exit  points  and  for  the 
specification  of  input  parameters  at  optional  levels  of  detail.  This  section  also  serves  as  a  reference 
which  explains  the  various  options  in  detail.  The  final  main  section  presents  the  card-image  input  data 
structure,  along  with  a  table  indicating  the  cards  required  by  the  different  options,  a  flowchart  relating 
input  and  output  data  with  the  subroutines  of  the  model,  and  a  sample  set  of  JSL  (Job  Specification 
Language)  commands  for  executing  the  program  on  the  ASC.  A  sample  execution  of  program 
COHORT  is  given  as  the  appendix. 
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HISTORICAL  SUMMARY  OF  PROPAGATION  IN  RANDOM  MEDIA 

The  modern  history  of  research  on  the  scattering  of  a  radiation  field  by  a  randomly  fluctuating 
continuum  began  during  and  immediately  after  World  War  II.  A  number  of  significant  studies  reported 
then  by  Bergmann  [6],  by  Mintzer  [7],  and  by  Pekeris  [8]  were  undertaken  to  explain  observations  of 
the  effects  of  temperature  fluctuations  in  the  ocean  on  a  propagating  acoustic  signal.  Since  the  tempera¬ 
ture  fluctuations  are  weak  and  the  range  of  the  experiments  of  interest  then  was  limited,  the  principal 
method  of  analysis  was  a  single-scatter  theory,  or  a  Born  approximation.  Subsequent  to  these  studies 
there  were  a  number  of  conceptually  and  mathematically  similar  ones  motivated  by  observations  of  the 
effects  of  atmospheric  turbulence  on  a  propagating  laser-beam  signal. 

Two  monographs,  by  Chernov  [9]  and  by  Tatarski  [10],  were  translated  into  English  and  pub¬ 
lished  in  1960  and  1961.  These  books  provided  a  rather  complete  survey  of  the  research  that  had  been 
carried  out  in  the  Soviet  Union,  and  they  defined  what  could  be  termed  the  state  of  the  art  at  that  time. 
The  research  discussed  pertained  to  both  acoustic  and  electromagnetic  radiation  and  emphasized  the 
importance  of  turbulence  as  the  dynamic  process  that  ultimately  gives  rise  to  the  scattering.  The 
manuscript  by  Tatarski  was  particularly  noteworthy  for  its  description  of  the  fluctuating  medium  and  for 
its  reliance  on  the  Kolmogorov  spectrum  as  a  correct  description  of  the  scattering  mechanism.  Both 
monographs  were  also  significant  for  introducing  the  Rytov  approximation  to  Western  researchers  of 
random-scattering  problems.  The  Rytov  approximation,  like  the  Born  approximation,  is  based  on  per¬ 
turbation  ideas,  but  the  claim  of  both  authors  was  that  the  approximation  correctly  accounted  for  the 
multiple-scattering  effects  necessary  for  the  calculations  to  be  valid  for  long  ranges. 

Research  into  the  subject  went  through  a  high  point  in  activity,  and  in  controversy,  during  the  six¬ 
ties,  with  most  of  the  reported  studies  treating  the  scattering  of  electromagnetic  signals  in  the  atmo¬ 
sphere.  Much  of  the  controversy  centered  around  two  questions:  What  measures  of  the  statistics  of 
the  radiation  field  are  most  conveniently  determined  in  physical  experiments  and  most  conveniently 
incorporated  in  theories?  How  does  one  derive  theories  that  properly  incorporate  multiple  scatter 
effects,  as  well  as  the  effects  of  diffraction?  A  number  of  studies  of  the  second  question  were  framed 
in  terms  of  the  relative  merits  of  the  Born  and  Rytov  approximations. 

By  the  close  of  the  decade  of  the  sixties  the  controversy  on  these  fundamental  questions  largely 
ceased.  The  central  role  of  the  multipoint  statistical  moments,  termed  coherence  functions  in  the  prop¬ 
agation  literature,  was  recognized  by  increasing  numbers  of  researchers;  techniques  had  been  developed 
for  deriving  theories,  in  the  form  of  differential  equations,  governing  these  statistical  moments;  specific 
equations  had  been  written  for  the  second-  and  fourth-order  moments,  the  most  crucial  moments  for 
discussing  experiments;  and  studies  were  frequently  reported  treating  the  analysis  and  the  solution  of 
these  equations  in  specific  applications.  A  second  monograph  by  Tatarski  (111,  which  appeared  in 
English  in  1971,  deemphasized  the  role  of  the  Rytov  approximation  highlighted  in  the  earlier  work 
The  position  espoused  in  this  second  monograph  appeared  to  be  quite  similar  to  that  reached  by  a  grow¬ 
ing  number  of  researchers  in  the  United  States. 

Research  in  the  seventies  addressed  the  need  to  solve  the  governing  field  equations  on  the 
second-  and  fourth-order  moments.  These  efforts  included  obtaining  analytic  solutions  for  idealized 
experiments  and  the  development  of  the  general  numerical  algorithms  needed  to  address  realistic 
experiments.  In  the  seventies,  also,  the  specific  application  appeared  to  shift  away  from  the  scattering 
of  laser  beams  in  the  atmosphere  to  the  scattering  of  acoustic  signals  in  the  ocean.  While  in  principle 
the  applications  are  the  same,  four  factors  make  an  ocean  acoustic  experiment  different  in  detail  First, 
the  specific  statistic  of  the  received  signal  of  interest  differs  for  some  acoustic  systems.  Thus,  for 
example,  estimation  of  the  temporal  spreading  of  an  acoustic  signal  received  by  a  single  omnidirectional 
hydrophone  is  of  interest  for  discussing  the  performance  of  communication  systems  for  which  the 
ocean  provides  the  channel.  This  can  be  distinguished  from  the  primary  motivation  for  much  of  the 
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work  on  electromagnetic  signals  carried  out  in  the  sixties:  the  derivation  of  algorithms  for  estimating 
the  performance  of  aperture  systems  for  which  the  most  important  measure  of  the  received  radiation  is 
of  the  directional  spreading  of  a  signal  across  the  face  of  an  optical  lens.  This  latter  application,  of 
course,  has  its  counterpart  for  signals  received  by  an  ocean  acoustic  array.  Second,  the  dynamic  process 
that  gives  rise  to  the  fluctuating  continuum  is  different  in  the  ocean  acoustic  application  [12].  Third, 
the  wavelengths  of  the  acoustic  signals  for  the  experiments  of  interest  are  large  (relative  to  important 
characteristic  lengths)  compared  to  previous,  nonacoustic,  applications.  Fourth,  the  ocean  is  an  inho¬ 
mogeneous  and  a  highly  anisotropic  propagation  medium  [13-15].  While  the  second  and  third  of  these 
four  factors  make  the  ocean  acoustic  experiment  different  from  the  previous  application  in  degree,  the 
first  and  fourth  can  make  it  different  in  kind.  It  is  clear  that  a  change  in  the  specific  statistic  that  is  of 
interest  would  change  the  nature  of  the  prediction  model.  It  is  less  clear,  although  equally  true,  that 
the  incorporation  of  inhomogeneity  and  anisotropy  of  the  medium  introduces  additional  length  scales 
that  need  to  be  parameterized  in  additional  nondimensional  ratios.  Valid  prediction  modeling  will  usu¬ 
ally  require  different  models  for  different  limiting  values  of  these  nondimensional  ratios. 

Continuing  with  a  synopsis  of  the  advances  made  in  the  seventies,  a  general  appreciation  of  the 
mathematical  identity  of  the  scattering  by  a  randomly  inhomogeneous  continuum  and  the  quantized 
motion  of  a  particle  in  a  randomly  perturbed  potential  field  [16,17],  as  well  as  a  duality  between  a  radia¬ 
tive  transport  theory  and  the  equation  governing  the  two-point  coherence  function,  was  achieved  during 
this  time  span  [18-20],  While  this  appreciation  has  not  greatly  altered  the  general  flow  of  the  develop¬ 
ment  of  theories,  it  has  introduced  new  techniques  that  could  prove  useful  in  solving  specific  problems, 
e.g.,  the  use  of  Monte  Carlo  calculations  [21].  Two  additional  analytic  techniques  were  introduced  into 
the  literature  of  stochastic  volume  scattering  in  the  seventies.  One  was  the  use  of  the  formalism  of 
Feynman  path  integrals  [22].  The  other  was  the  formulation  of  the  scattering  problem  in  terms  of  a 
modal  expansion;  the  scattering  mechanism,  in  this  formulation,  results  in  a  coupling  via  intermodal 
energy  transfer  of  the  normal  modes  defined  for  a  depth-dependent  background  medium  (23-26).  The 
motivation  for  the  modal  expansion  formulation  was,  clearly,  a  realization  that  the  ocean  environment 
does  define  a  waveguide,  which  becomes  more  obvious  with  the  decreasing  frequencies  dictated  by 
changing  applications. 

A  number  of  extended-study  programs  of  the  random  scattering  of  acoustic  signals  by  tempera¬ 
ture  fluctuations  were  carried  out  throughout  the  seventies.  Perhaps  the  most  extensive,  and  certainly 
the  most  extensively  reported,  were  the  studies  of  the  JASON  group;  a  readable  summary  of  much  of 
their  effort  has  been  published  in  book  form  [27],  There  appear  to  have  been  four  principal  contribu¬ 
tions  of  the  JASON  studies.  First  was  the  emohasis  they  placed  on  the  need  to  relate  the  acoustic 
event,  the  stochastic  scattering,  to  the  oceanographic  events,  the  presence  of  interna)  waves  as  the  con¬ 
trolling  dynamic  process  and  of  a  depth -dependent  background  sound-speed  profile.  Second  was  the 
importance  of  combining  important  characteristic  length  scales  into  nondimensional  parameters,  which 
can  then  be  used  to  classify  scattering  experiments  according  to  separate  domains  of  parameter  space. 
Since  the  JASON -group  research  accepts  the  validity  of  a  model  that  is  rigorously  derived  for  a  medium 
that  is  both  homogeneous  and  isotropic,  only  two  nondimensional  parameters  are  required  to  classify  all 
scattering  experiments:  one  is,  essentially,  the  ratio  of  the  experiment  range  to  the  range  at  which 
diffraction  effects  become  significant,  and  the  other  is,  essentially,  the  ratio  of  the  experiment  range  to 
the  range  at  which  significant  acoustic  energy  (say  one-half  of  the  original  energy)  has  been  scattered. 
The  third  contribution  was  the  introduction  of  the  Feynman  path  formalism  already  alluded  to  above. 
The  fourth  contribution  was  a  detailed  discussion  of  a  number  of  reported  experiments,  principally 
those  reported  by  Ellinlhorpc  et  al.  (28)  and  by  Ewart  [29).  Although  the  basic  formalism  presented  by 
the  JASON  group  can  be  applied  to  a  broad  spectrum  of  experiments,  it  has  been  applied  in  detail  only 
to  experiments  in  which  the  spatial  resolution  of  the  signal  was  not  a  principal  objective.  All  of  the 
comparisons  considered  only  the  statistics  of  a  time  series  measured  at  a  single  point  in  the  acoustic 
field. 
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A  second  extended  program  of  studies  was  carried  out  by  a  less-well-defined  group,  centered  pri¬ 
marily  either  at  New  York  University  or  around  Tapper!  (see,  for  example,  Refs.  16-18,  21,  25,  and  26). 
The  scope  of  the  program  carried  out  by  Tappert  et  al.  was  more  limited  than  that  of  the  JASON  group; 
Tappert  was  essentially  interested  in  the  lower  order  spatial  statistics,  oi  in  the  "aperture  problem"  that 
motivated  much  of  the  electromagnetic  work  of  the  fifties  and  sixties.  Further,  this  second  effort  was 
exclusively  either  analytical  or  numerical  in  nature;  no  reference  was  made  to  any  specific  series  of 
experiments.  The  principal  contributions  were  noted  earlier  in  this  synopsis:  namely,  elucidation  of  the 
relationship  between  the  ocean  acoustics  problem  and  that  of  the  quantized  motion  of  a  particle  in  a 
perturbed  potential  field  116,1  /J,  the  use  of  Monte  Carlo  calculations  (21],  and  the  formulation  of  the 
scattering  problem  in  terms  of  a  modal  expansion  (25,26). 

The  third  program  of  studies  was  carried  out  by  researchers  at  NRL  and  by  Beran  and  McCoy, 
working  with  NRL.  Once  again  the  scope  of  the  study  program  was  limited  to  estimation  of  the  spatial 
coherence  across  a  receiving  aperture.  There  were  four  principal  accomplishments.  The  first  was  the 
demonstration  that  the  degree  of  anisotropy  of  the  scattering  mechanism,  i.e.,  that  due  to  internal 
waves,  was  such  as  to  necessitate  the  introduction  of  a  new  nondimensional  ratio  for  a  complete  param¬ 
eterization  of  ocean  acoustic  experiments  (13-15).  The  scattering  models  presented  in  the  optical  litera¬ 
ture  of  the  sixties,  or  by  the  JASON  group  in  the  seventies,  can  be  shown  to  be  valid  as  a  limit  for 
small  values  of  this  anisotropy  ratio.  This  limit  can  be  justified  either  for  propagation  in  an  isotropic 
medium  or  for  the  propagation  of  high-frequency  signals.  (For  typical  experiments,  high  frequency 
implies  greater  than  several  hundred  hertz.) 

Second,  the  initial  analysis  by  the  NRL  group  was  concerned  with  making  a  number  of  approxi¬ 
mations  in  order  to  develop  a  closed-form  expression  for  estimating  the  loss  of  spatial  coherence  over  a 
billboard  array  (13,30).  The  purpose  of  this  development,  which  forms  the  basis  of  the  computational 
algorithm  in  the  present  report,  was  to  provide  a  realizable  calculation  which  could  be  programmed  for 
antenna-performance  estimates  of  ail  types.  The  many  assumptions  requited  were  to  be  physically 
based  and  tested  with  experiments  in  the  ocean  This  effort  has,  to  the  best  of  our  knowledge,  turned 
out  to  be  successful  for  horizontal  antennas,  in  that  an  extensive  series  of  experimental  results  compar¬ 
ing  favorably  with  the  estimates  of  horizontal  coherence  produced  with  the  program  described  here  has 
been  assembled  and  is  being  published  The  initial  effort  for  combined  horizontal  and  vertical  apertures 
showed  that  multiple  parameters  were  required  Further,  it  was  judged  most  practical  to  produce 
revised,  simplified,  and  improved  versions  with  more  limited  scope:  separate  horizontal  and  vertical 
antenna  models,  for  example  l  he  first  such  circumscribed  effort  addressed  the  vertical  coherence  or 
diffusion  ror  a  single  path.  T'os  effort  was  carried  to  conclusion  and  is  available  in  a  calculalional  pack¬ 
age,  COVKRT.  which  produces  estimates  of  single  path  scattering  along  the  selected  macroray  (1)  This 
calculation  is  for  those  circumstances  wherein  appeal  to  the  simplification  of  a  discrete  field  characteri¬ 
zation.  rays,  is  warranted  These  cases  are  basically  short  range,  where  the  multipath  effects  can  be 
clearly  separated  from  tnc  spreading  in  the  vertical  angle  estimated  by  this  program. 

The  third  contribution  of  the  NRL  group  has  been  the  collection  of  a  large  body  of  experimental 
material  renting  tc-  the  diffusion  of  an  acoustic  field  in  the  norizontal  regime,  as  measured  for  horizon¬ 
tal  line  arrays,  and  the  comparison  of  these  data  with  predictions  from  the  closed-form  solutions.  Dur¬ 
ing  the  formulation  phases  of  these  efforts,  many  lengthy  discussions  were  held  regarding  the  manner 
in  which  component  coherences  of  a  multipath  structure  combine  to  form  the  coherence  of  the  total 
field  For  horizontal  antennas,  however,  all  the  results  to  date  indicate  that  the  incoherent  summation 
of  the  coherences  of  multipath  components  ts  a  feasible  simplification  A  similar,  although  esscr.  body 
of  experimental  material  has  been  assembled  and  published  concerning  the  distribution  of  the  vertical 
multipath  field 

Ti  c  fourth  stage  of  the  NRL  group’s  analysis  has  concerned  a  combined  estimation  of  the  coher¬ 
ence  effects  of  refraction  and  scattering,  the  model  used  is  called  the  Combined  Kffccts  Model  (CKM) 
(31  In  this  calculation  the  simplifications  inherent  in  ihc  assumption  of  a  sir.^le-paih  field  arc  removed 
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and  the  coherence  function  itself,  with  an  implicit  smoothed  multipath  field,  is  propagated  through  the 
medium.  This  calculation  combines  deterministic  (multipath)  and  stochastic  (scattering)  coherence 
mechanisms,  and  it  gives  the  solution  to  vertical  multipath  problems  that  are  calculated  to  ranges  which 
render  the  discrete-field  model  of  COVERT  unacceptable.  Whether  the  discrete-field  vertical  model  or 
the  more  complicated  vertical  model  termed  CEM  is  appropriate  for  a  given  situation  depends  upon  the 
objective  of  the  study. 

In  summary,  low-frequency,  transverse  acoustic  dispersion  in  the  deep  ocean  has  been  measured 
and  modeled  by  the  NRL  group.  Transverse  horizontal  coherence  resulting  from  stochastic  volume 
inhomogencity  has  been  most  intensively  studied,  and  it  is  the  subject  of  the  present  prediction-model 
reported,  COHORT.  Transverse  vertical  coherence  has  been  treated  in  two  ways,  as  described,  and  will 
be  covered  in  subsequent  reports. 

ANALYTIC  BACKGROUND 

In  the  Introduction,  several  concepts  basic  to  understanding  the  COHORT  program  were  intro¬ 
duced.  In  this  section,  we  define  these  basic  concepts  and  show  the  strengths  and  limitations  by  means 
of  examples. 

First,  the  nature  of  the  oceanic  environmental  anistropy  and  '.he  source-receiver  locations  clearly 
defines  three  distinguishable  directions,  as  illustrated  in  Fig.  1:  the  depth  direction  r,  a  horizontal  direc¬ 
tion  determined  by  the  source/ receiver  locations  r,  and  a  horizontal  direction  transverse  to  that  deter¬ 
mined  by  the  source/receiver  locations  y. 


Eif.  I  —  CoorUuutc  tyucra 


The  mathematical  measure  of  the  eo  rence  of  an  acoustic  field,  which  is  the  fundamental  con¬ 
cept  of  this  report,  is  the  ensemble-averaged  product  of  the  complex  pressure  field  and  its  conjugate 
measured  at  two  points.  Wc  write 

<!'(«,.  x,)>  -  <p(x,)p»(x j)>.  (I) 

where  <F(x,  ij)>  denotes  the  coherence  function,  an  asterisk  denotes  cc  tplex  conjugation,  and  the 
angular  brackets  denote  ensemble  a  vet  aging. 

Gur  first  example  is  designed  to  show  that,  although  the  coherence  function  is  a  stochastic  formu- 
la.ion,  >;  is  influenced  by  purely  deterministic  environmental  variations,  as  well  as  by  random  ones.  Wc 
consider  an  acoustic  field  given  by  two  plane  waves,  each  propagating  in  the  x-z  plane,  independent  of 
cross  range  y,  in  directions  that  arc  inclined  to  the  x  axis  by  angles  d>(  and  <Jj,  respectively.  (The 
geometry  is  shown  in  Fig.  2.)  An  experiment  in  which  the  received  signal  can  be  resolved  into  arrivals 
from  two  different  paths  of  a  multipath  structure  is  approximated  by  this  calculation.  The  complex 
pressure  field  for  the  two  plane  waves  is  written 
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Fig.  2  —  Angles  of  vertical  multipath 


2 


p(x)  ■»  Ax  exp  (dc(sin  d>xz  +  cos  d^x)]  +  A 2  exp  t/Jc  (sin  4>2z  +  cos  4>2x) ),  (2) 

where  k  -  2ir/\  -  w/c  is  the  signal  wavenumber,  c  being  the  sound  speed  in  the  medium.  We  form 
<T(i  +  s/2,  z  -  s/2,  x)>  for  two  points  (z  ±  s/2),  separated  by  a  distance  s,  in  the  same  range 
plane  (x).  for  an  experiment  in  which  A\  and  A2  are  deterministic  (complex)  constants.  The  result  can 
be  written  as  a  sum  of  three  terms,  i.e., 

<r(z  +  s/2,  z  -  s/2.  x)>  -  <r,(z.  s.  x)>  +  <r2(z.  s.  x)>  +  <rtu  s.  x)>,  o) 

where 

<r,(z.  s.  x)>  -  |d„iJcxp  (ik  sin  it  -  l.  2.  (4a> 

and 


<I\(z.  s.  x)>  -  2  Re  |-4t/4j*cxp  (cfc [(sin  -  sin  </>2)z 

,  , ,  i  [(  sin  <^t  +  sin 
+  tcos  d>i  -  cos  6>)x))j  exp  I/Jcj - - - 


.  (4b) 


The  functions  <Tj>  arc  the  components  of  the  totai  <I‘>  that  can  be  associated  with  the 

two  plane  waves  existing  independently  of  one  another;  the  function  <  l't>  is  an  interference  term. 


The  ‘direct  coherence*  functions  (<!')>  and  <!'»>)  oscillate  with  changes  in  the  separation 
coordinate  (i),  and  these  osculations  arc  a  result  of  the  teceiver  pc >.ts‘  lying  on  a  plane  that  makes  an 
angle  with  the  phase  planes  of  tht  incoming  waves.  It  is  well  appreciated  that  the  osc*Hation«  of  either 
one  of  the  direct  coherence  functions  arc  readily  removed  by  introducing  a  phase  shift  in  the  signal 
received  at  one  of  the  points,  that  is.  by  steering  the  array  to  the  direction  of  the  incoming  signal. 


The  ‘mteiference  coherence’  function  (<!%>)  oscillates  with  changes  in  either  the  coordinates 
<x.  i)  or  the  separation  coordinate  (if  Further.  it  is  clear  from  the  defining  equation  the. t  these  oscilla¬ 
tions  cannot  hi  removed  by  the  introduction  of  a  phase  shift  in  the  stgna!  received  at  one  of  the  end 
points  (It  requires  an  averaging  over  all  possible  pair  separations  of  a  ‘filled’  array  to  remove  *  , *  >  . 
see  below  )  Oscillations  in  pairwise  signal  coherence  with  changes  «n  the  positions  of  the  receiver 
points  must  always  be  expected  in  an  experiment  <n  which  the  received  ogne.  arrives  from  more  than 
one  direction  They  are  a  manifestation  of  a  deicrmir.Mtc  phenomenon,  i  e  .  interference,  and  are  not  to 
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be  confused  with  a  "loss"  of  signal  coherence  that  results  from  statistical  fluctuations  in  a  single-path 
field. 


If  one  transforms  the  coherence  function  with  respect  to  separation  (s)  to  obtain  a  quantity 
measurable  by  an  array,  we  find  that  the  interference  term  <TC>  gives  rise  to  an  expression  that  oscil¬ 
lates  rapidly  as  a  function  of  the  depth.  The  following  statements  can  be  demonstrated: 

•  For  an  array  sufficiently  long  to  resolve  the  directions  of  the  two  incoming  plane  waves,  the 
smoothing  or  averaging  of  values  over  the  z  coordinate  will  be  sufficient  to  cancel  the  oscillations. 

•  For  an  array  that  is  insufficiently  long  to  resolve  the  directions  of  the  two  incoming  plane 
waves,  the  averaging  of  values  will  not  be  sufficient  to  cancel  the  oscillations.  The  nonzero  averaged 
value  of  the  expression  in  this  case  is  a  measure  of  the  interference  of  the  unresolved  plane  waves. 


The  main  point  of  the  preceding  summary  was  to  emphasize  that  there  are  a  number  of  "deter¬ 
ministic"  factors  that  affect  the  pairwise  coherence  of  a  narrowband  acoustic  field  as  a  function  of 
separation  distance  and  that  these  factors  must  be  considered  before  any  others.  Among  these  deter¬ 
ministic  factors  are  the  orientation  of  the  measuring-point  pair  set  to  a  single  propagation  direction,  as 
well  as  the  number  of  components  and  the  angular  separation  of  a  multipath  structure.  We  refer  to 
issues  of  this  type  as  deterministic  coherence  effects. 


Randomness  can  enter  discussions  of  coherence  in  two  ways.  One  way  is  to  introduce  random 
perturbations  of  the  parameters  that  describe  the  deterministic  factors  discussed  above.  These  random 
effects  on  coherence  are  usually  overridden  by  the  geometric  and  algebraic  aspects  of  the  definition  and 
will  not  be  considered  further.  A  more  inclusive  aspect  of  randomness  can  be  introduced  by  allowing 
the  phase  functions  of  the  complex  amplitudes  (Ax  and  A2)  to  be  centered  stochastic  processes  defined 
over  the  spatial  coordinate.  Thus,  the  expression  for  p(x)  given  by  Eq.  (2)  is  replaced  by 

p(x)  -  Mil  exp  [ik(s\n  $xz  +  cos  4>ix)  +  i<t>i(z,  x)l 

+  M2|  exp  {/A (sin  <£2z  +  cos  4>i x)  +  /<J>2(z,  x)l,  (5) 


where  d>i  and  <l>2  are  residual  stochastic  phase  functions.  We  form  <T(z  +  s/2,  z  -  s/2,  x)>  as 
before,  and  the  result  can  again  be  written  as  the  sum  of  three  terms,  as  in  Eq.  (3),  where  now 

<r„(z,  s,  x)>  -  \A„\2  exp  ( ik  sin  <exp  / [4»„ (z  +  s/2,  x)  -  <Mz  -  s/2,  x)]>  (6a) 


and 


<F c(z,s,x)>  -  U,|  U2|  exp  ik  — "  +  S'n 


2 


x  (  B  exp  i 


<t>. 


z  + 


*2 


S 

z-  2 


+  B*  exp  — / 


z  - 


<t>2 


z  + 


(6b) 


where 


B  =  exp  { /A {(sin  |  -  sin  <j>2)z  +  (cos  4>\  -  cos  </>2)x]}.  (6c) 

To  simplify  the  averaging  in  Eq.  (6)  it  is  necessary  to  say  something  of  the  statistics  that  govern  the  <!>„. 
For  illustrative  purposes  only  we  assume  homogeneous,  Gaussian  statistics,  an  assumption  that  enables 
us  to  write 

<r*(z,  s.  x)>  =  M„|2  exp  (ik  sin  <f>ns)  exp  {— l/2[o- 2(0)  -  <r2(s)]J,  n  -  1,  2,  (7a) 


8 


NRL  REPORT  8oi5 


and 


<rr(z,  s,  x)>  **  2  R.e  |  [  /4 1  /4  2 I  exp  { /A:  I  (sin  <t>\  -  jin  <t>2)z  +  (cos  <f>\  -  cos  fo)-*]]] 

exp  [  -  1/2  (a2(0)  +  a2(0)  —  2o- j22(5) j) , 


x  exp 


sin  <f>i  +  sin  <t>2 
ik  - - - s| 


where 


cr j(s)  =  <4>n(0)<I>„(s)>,  n  =  1,  2, 
are  the  correlation  functions  of  the  two  ptocesses,  taken  separately,  and 

o-|J2(s)  =  <<t>1(0)4>j(s)> 


(7b) 

(8a) 

(8b) 


is  the  cross-correlation  function  of  the  two  processes,  taken  jointly.  For  physically  realizable  systems, 
the  values  of  the  correlation  and  cross-correlation  functions  are  greatest  for  zero  separation  distance  end 
are  zero  for  infinite  correlation  distance.  Also,  2a22(0)  <  cr?( 0)  +  a2  ( 0)  for  physically  realizable  sys¬ 
tems. 


Comparing  Eqs.  (7)  to  Eqs.  (4),  we  see  that  the  introduction  of  a  randomly  varying  phase  results 
in  an  exponential-decay  term,  with  an  exponent  that  depends  on  the  separation.  It  is  this  decay  term  that 
describes  a  coherence  loss  that  is  an  inherent  feature  of  the  statistical  nature  of  the  acoustic  field  and  is 
representative  of  the  stochastic  coherence  effect  of  interest  in  this  report. 


Wr  *  ».*  --  >er  Representations  of  Acoustic  Fields 

The  acoustic  coherence  function  has  been  deSned  by  Eq.  (1).  For  the  remainder  of  this  report  we 
will  consider  horizontal  coherence  only.  That  is,  *i  and  x2  will  be  restricted  to  a  horizontal  line  at  a 
depth  z0  and  a  range  x.  This  line  will  define  the  y-axis.  The  acoustic  field  is  assumed  to  be  homogene¬ 
ous  in  the  (horizontal)  >*-direction,  in  which  case  <F>  is  independent  of  y.  Therefore,  the  two-point 
transverse-horizontal  spatial  coherence  can  be  written  <r(s/,  x,  z^>.  In  the  following,  it  will  be  writ¬ 
ten  simply  as  r(s),  with  the  dependency  on  x  and  z0  and  the  ensemble  averaging  assumed  to  be  under¬ 
stood. 

Analogous  to  the  time/frequency  Fourier  transform  pair  is  the  one-dimeasional-spaee/ wave- 
number  transform  pair, 

g(k)  -  f  ^  f^e'^dy  (9a) 

and 

fly)  -  glk)e*'^-^-.  (9b> 

2tr 

If  we  define  F  to  be  the  two-point  coherence  of  the  signal-only  portion  of  the  acoustic  field,  the 
wavenumber  spectrum  of  the  signal  (at  a  fixed  frequency  a>0)  is 

Slk )  -  f  ’  -k0  k  ^  *0.  (10) 

-a* 

where  k0  «■  t»Jc  and  c  is  sound  speed. 

The  spectrum  Slk)  represents  the  (horizontal)  arrival  structure  of  the  signal.  This  is  seen 
through  the  relationship 

k  -  k0  sin  6.  (11) 
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A  horizontal  receiving  line  array  is  assumed  to  lie  along  the  .p-axis,  in  which  case  the  angle  9  is  mea¬ 
sured  relative  to  broadside  of  the  array  (see  Fig.  1).  For  example,  a  perfectly  coherent  plane  wave  of 
unit  amplitude  arriving  broadside  to  a  linear  array  has  a  spatial  coherence  function  (along  the  array)  of 

T(s)  =  l  (12a) 


and  a  wavenumber  spectrum  of 

S(k)  =  S(k0  sin  9)  =  5(0)  •  2 it,  (12b) 

where  S(-)  is  a  Dirac-delta  function.  If  we  drop  the  constant  ka  and  consider  the  spectrum  as  a  func¬ 
tion  of  sin  6,  the  spectrum  is  zero  except  when  sin  9  =  0  (i.e.,  B  =■  0,  which  corresponds  to  broadside). 
Figure  3  shows  a  plot  of  Eq.  (12b).  For  a  perfectly  coherent  unit-amplitude  plane  wave  arriving  from 
30°  off  broadside. 

r(s)  -  e'*otoW6>s  (i3a) 

and 

S(k)  =  2rr8  (k  —  k0  sin  ir/6)  -  2n8  k0( sin  6  -  sin  it/ 6).  (13b) 


Figure  4  shows  a  plot  of  Eq.  (13b). 


T 

-i 


0 

Fig.  3  —  Plot  of  Eq.  (12b) 


— | - sinfl  . 

+  1 


-1 


T  sin7r/6  1 
0  +1 


sinfl  . 


Fig.  4  —  Plot  of  Eq.  03b) 


The  wavenumber  spectrum  of  a  partially  coherent  signal  field  is  not  a  delta  function,  but  rather 
the  spectrum  is  spread  in  angle,  as  shown  in  Fig.  5,  Thus,  the  spectrum  S(k)  is  seen  to  represent  the 
(horizontal)  arrival  structure  of  the  signal  as  a  function  of  arrival  angle  9.  Furthermore,  8  varies 
between  -ir/2  and  tt/2  from  broadside  (and  simultaneously  between  rr/2  and  3tt/2,  the  front/back 
ambiguity).  Also,  for  propagating  waves  (as  opposed  to,  say,  electronic  noise), 

S(k)  -  0  for  U I  >  kq.  (14a) 


or  equivalently. 


S(k)  -  0  for  (sin  fll  >  1. 


(14b) 


in  which  case,  the  transform  pair  may  be  written 

S(k)  —  f  r(s)c"‘*‘ds 


and 


(!5u) 

(15b) 
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sin(J 


Fig.  5  —  Wavenumber  spectrum  of  partially  coherent  signal  field 


Transverse  Horizontal  Coherence  Model 

A  formulation  for  the  transverse  spatial  coherence  of  an  acoustic  wavefront  in  an  inhomogeneous 
and  anistropic  ocean  has  been  developed  under  the  auspices  of  NRL  [13-15,19,20,30].  A  reduced  form 
estimates  horizontal  coherence  generated  by  three-dimensional  spatial  departures  of  temperature  from 
its  mean  and,  this,  is  directly  suitable  for  long,  horizontal  line  arrays.  The  reduced  equation  for  nar- 
ro'jban<_  signal  coherence  along  a  horizontal  line  transverse  to  the  direction  of  (entirely-waterborne) 
propagation  from  a  point  source  is  [14] 

f(s)  -  F(0)  exp  [-£fr0i/2xs3/2],  (16) 


where 

x  =  range  from  the  source  to  the  receiver  (m), 
k0  =  acoustic  wavenumber  (m-1),  and 
s  =  horizontal  separation  of  measurement  points  (m); 


and 


where 


£=  1.7 


IvAj, 


average  sound  speeu  alo"g  path  of  signal  propagation  (m/s), 
temperature  (°C), 

average  derivative  of  sound  speed  with  respect  to 
temperature  along  path  of  signal  propagation  f(m/s)/0C], 


ly  «*  vertical  correlation  leng'h  of  tempera'u.e  fluctuation*  (m),  and 
A}  =  nominal  strength  of  random  temperature  field  (°C2/m). 

The  above  parameters  are  discussed  lusher  in  a  later  sect.on. 


The  normalized  coherence  function  is  defined  by 

C(s)  =  r(s)/r(o) 

=  exp  {—  1konxsin] . 


(17) 


The  corresponding  wavenumber  spectrum  of  the  signal  is 

5(A)  =  r(0)  exp  (-£A05/2xs3/j)  exp  (-iks)ds  (18) 

=  r(0)5(A). 

Generally,  the  spectrum  S(k)  is  determined  by  numerical  integration  (transform)  because  of  the  3/2 
exponent  on  the  spatial  variable. 
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The  range  of  validity  for  the  transverse  coherence  model  was  estimated  [30,  p.  14]  to  be  100  Hz  to 
300  Hz.  A  recently  published  set  of  data,  however,  shows  that  it  performs  well  from  100  Hz  to  400  Hz. 
Further,  the  model  estimates  were  developed  as  harmonic-wave  solutions,  i.e.,  narrowband  solutions. 
However,  since  dispersion  is  estimated  to  be  minor,  the  estimates  are  expected  to  be  valid  for  any  prac¬ 
tical  low-frequency  signal  bandwidth  [31,  Fig.  5], 


Beam  Patterns 


Similarly  to  the  wavenumber  spectrum  of  the  acoustic  field,  the  response,  or  array  pattern,  of  a 
horizontal  receiving  line  array  has  a  wavenumber  representation.  The  array  pattern  (at  fixed  frequency 
coq)  of  a  linear  array  of  n  equally  spaced,  omnidirectional  elements  is  given  by 

^  12 


B0(&) 


sin  ( nk0  —  sin  9) 


n  sin  (&0  —  sine) 


(19a) 


where  d  is  the  element  spacing.  In  wavenumber  space, 

sin  ( nkd/2 ) 


B0(k) 


n  sin  Vcd/2) 


(19b) 


The  above  equations  hold  only  for  equally  spaced,  perfectly  linear  (nondeformed)  arrays.  If  the  ele¬ 
ments  are  not  omnidirectional  but  have  identical  beam  patterns  Bs(9),  the  total  (or  product)  array  pat¬ 
tern  is 


B„(k)  -  Ba(k)Bg(k), 

where  Ba(k)  is  the  array  pattern  that  would  result  if  all  the  elements  were  omnidirectional.  A  practical 
example  is  one  in  which  each  identical  element  of  the  array  consists  of  a  group  of  omnidirectional 
hydrophones  equally  spaced  along  the  array  axis.  In  this  case,  the  element  pattern  Bg(k)  takes  on  the 
form  of  B0(k ),  where  now  n  is  the  number  of  omnidirectional  hydrophones  in  each  group,  and  d  is  the 
spacing  between  individual  group  hydrophones.  Now, 

B„{k )  =  Ba(k)B0(k). 


Normally,  the  outputs  from  the  different  hydrophones  of  a  group  are  summed  (e.g.,  in  the  array) 
without  time  delays.  That  is,  the  element  or  group  pattern  Bg(k)  is  permanently  steered  to  broadside. 
The  receiving  pattern  is  steered  by  time-delay  or  phase-rotation  of  the  outputs  from  the  different  group 
elements  (i.e.,  by  the  steering  of  the  array  pattern).  For  example,  if  the  array  pattern  is  steered  in  the 
direction  k ,  *=  sin  6  ,  the  product  pattern  is 

Bp (k,  kr )  =  Ba(k-  kr)  Bg(k). 

Note  that  the  mqjor  response  of  this  pattern  is  not  necessarily  in  the  direction  kr. 


Array  Signal  Gain 

For  a  signal  arriving  from  direction  9S  ( ks  «  ko  sin  9S),  with  the  receiving  line  array  steered  to  the 
direction  9„  the  power  out  of  the  beamformer  due  to  the  signal  only  is  given  by 

P(kr,  ks)  -  (]k°  Bp^k’  *r>  Slh~  *.)  (21) 

%/~'Ko  Ztt 

-  r(0)  /  Ba{k  -  kr)  Bs  (. k )  S  (k  —  ks) 
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For  a  perfectly  coherent  plane  wave, 

S(k)  =  S (*  -  ks)  •  lit 
and 

PQ(kr,  ks)  =  r(0)fia(/cs  -  kr)  Bg(ks).  (22a) 

When  the  array  is  steered  to  a  signal  arriving  broadside, 

/yo,  o)  =  r(o)Bfl  (o)^(o).  (22b) 

Assuming  that  the  group  pattern  is  well  known,  it  can  be  normalized  such  that 

”  1-  (23a) 

Furthermore,  the  array  pattern  can  be  normalized  such  that  under  design  conditions  (e.g.,  the  array  is 
perfectly  linear) 


Ba(0)~l  (23b) 

[otherwise,  Ba(0)  <  1],  In  this  case, 

P0( 0,  0)  -  T(0).  (23c) 

When  the  array  is  steered  to  a  signal  arriving  from  the  off-broadside  direction  ks , 

P0(ks,  ks)  -  r(Q)Ba(0)Bg(ks)  (24) 

-  P0( 0.  0 )Bg(ks). 

Degradation  in  array  signal  gain  is  defined  by 

L(kr,  k,)  ~  P(k„  ks)/P0(Q,  0),  (25) 

where  P0(0,  0)  is  evaluated  under  design  conditions  [i.e,  fla(0)  -  1,  Fo(0,0)  =  r(0)].  Therefore,  for 
an  arbitrary  signal  spectrum, 

L(k„  ks)  =  J*  *°  Ba(k  -  k,)  Bg  (k)  S  (k-  ks)  (26) 

0  ITT 

When  the  array  is  steered  to  a  perfectly  coherent  plane  wave  arriving  broadside, 

L0(0.  0)  -  fla(0).  (27) 

This  represents  the  loss  (if  any)  due  to  the  main-axis  response  of  the  array  pattern  (due,  say,  to  the 
array  becoming  nonlinear  and  Ba(0)  <  1). 


If  the  array  is  steered  to  a  perfectly  coherent  plane  wave  arriving  off  broadside, 

Lo(k„  k,)  -  Ba(0)Bg(ks),  (28) 

and  an  additional  loss  factor  of  Bg(k,)  is  incurred  (due  to  the  group  pattern).  The  loss  factor  Ba(0) 
cannot  be  corrected  for,  since  it  represents  a  loss  due  to  a  deviation  from  expected  conditions 
[2?a(0)  *  1],  However,  for  any  steering  direction  kn  Bg(kr)  is  assumed  known  and  hence  can  be 
corrected  for  by  scaling  of  the  product  pattern,  i.e.. 


Bp(k.  k,)  -  Ba(k  -  k,)  Bg  (k)/B'(k,) 


which  ensures  that 

ks)  -  Ba( 0) 

when  the  array  is  steered  to  direction  ks.  Accordingly,  wc  rewrite  Eq.  (26)  as 


L(kr, 


1 


BAkr) 


S- 


°  Ba(k  -  kr)BAk)S(k  -  k.)  —■ 

o  2n 


(29a) 

(29b) 

(30) 
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MODEL  DESCRIPTION 


The  COHORT  model  is  an  implementation  of  the  theory  presented  in  the  previous  section.  Its 
organization  consists  of  two  basic  parts.  The  first  part  is  concerned  with  computing  signal  coherence 
and  the  second  with  subsequent  array  performance.  Furthermore,  there  are  three  primary  options  for 
how  the  COHORT  model  may  be  executed.  These  options  are  the  following: 

®  Compute  signal  coherence  only. 

•  Compute  signal  coherence  and  subsequent  array  performance. 

•  Input  previously  calculated  signal  coherence  data  and  compute  array  performance. 

The  flexibility  of  the  model  also  allows  for  program  entry  and  termination  at  different  stages  within 
these  options  and  allows  input  data  to  be  specified  in  varying  degrees  of  detail. 

COHORT  consists  of  a  driver  program  which  calls  a  sequence  of  subroutines  in  accordance  with  a 
user-defined  task.  Briefly,  these  subroutines  are  as  follows: 


MEDIUM  Inputs  and/or  computes  environmental  parameters  needed  to  evaluate  the  spatial- 
coherence  model. 

SPACOH  Computes,  prints,  and,  optionally,  plots  and  writes  the  signal  spatial-coherence  function 
onto  an  output  file. 

WAVCOH  Computes  and  prints  the  signal-coherence  spectrum  of  either  a  computed  or  an  input 
spatial-coherence  function.  The  spectrum  is,  optionally,  plotted  and/or  written  onto  an 
output  file. 

BEAMS  Computes  and,  optionally,  plots  the  group  pattern.  Reads  in  or  computes  and,  option¬ 
ally,  plots  the  array  pattern.  Reads  in  the  number  and  amounts  of  array  steers  or  shifts 
(i.e.,  signal-arrival  angles)  and  sets  up  information  for  convolution. 

CONVLV  For  each  array  pattern  shift  (signal  angle),  the  product  array  pattern  is  computed  and 
convolved  with  the  signal-coherence  spectrum  to  calculate  losses  in  array  signal  gain, 
which  are  printed  and,  optionally,  plotted. 


Figure  6  indicates  the  possible  entry  and  exit  points  in  the  above  sequence  of  subroutines. 


r 


SIGNAL  COHERENCE 


r 


ARRAY  PERFORMANCE 


Fig.  6  -  Entry  and  exit  options  of  program  COHORT 
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Program  COHORT 


Program  COHORT  is  a  driver  program  which  calls  a  sequence  of  subroutines  in  accordance  with 
task  specified  by  the  input  variables  KBEAM,  LCOH,  and  KCOH.  The  parameters  LCOH  and  KCOH 
refer  to  input  and  output  file  unit  numbers,  whereas  KBEAM  determines  the  primary  function  to  be 
performed.  There  are  three  cases  definable  by  the  value  of  KBEAM,  as  follows: 

Case  1 

KBEAM  -  1:  The  spatial  coheren  c  model  is  evaluated  and  the  coherence  spectrum  of  the 
signal  is  computed.  Then  the  resulting  Jesses  in  array  signal  gain  are  determined. 

If  LCOH  >  0,  the  spatial-coherence  function  is  written  onto  unit  LCOH. 

If  KCOH  >  0,  the  coherence  spectrum  is  written  onto  unit  KCOH. 

Case  2 

KBEAM  -  2:  Either  the  spatial  coherence  function  is  read  in  and  its  spectrum  computed,  or 
the  coherence  spectrum  is  read  in.  In  either  case,  losses  in  signal  gain  are  computed. 

If  LCOH  >  0,  the  spatial-coherence  function  is  read  from  unit  LCOH  and  its  spectrum  com¬ 
puted.  If,  also,  KCOH  >  0,  the  spectrum  is  written  onto  unit  KCOH. 

If  LCOH  <  0,  the  coherence  spectrum  is  read  from  unit  KCOH  (must  be  positive). 

Case  3 


KBEAM  -  3:  No  losses  in  signal  gain  are  computed.  The  NRL  coherence  model  is  evaluated 
and  the  computation  of  its  spectrum  is  optional. 

If  LCOH  >  0,  the  evaluated  spatial-coherence  function  is  written  onto  unit  LCOH. 

If  KCOH  >  0,  the  coherence  spectrum  is  computed  and  written  onto  unit  KCOH. 

If  KCOH  —  0,  the  coherence  spectrum  is  not  computed. 

If  KCOH  <  0,  the  coherence  spectrum  is  computed  but  not  written  to  an  output  file. 

Signal-Coherence  Subroutines 

There  are  three  subroutines  concerned  with  computing  signal  coherence:  MEDIUM,  SPACOH, 
and  WAVCOH. 

Subroutine  MEDIUM 

This  subroutine  determines  the  parameters  E,  k.  and  x  of  the  NRL  coherence  model.  Subroutine 
MEDIUM  is  activated  only  if  the  evaluation  of  the  signal  spatial-coherence  function  is  requested.  In 
subroutine  MEDIUM,  the  variables 

R,  VL,  CV,  DRT,  A,  and  E 
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are  read  in  a  card-image  format,  where 


R  =  x 
VL  =  lv 
CV=  c 

tie 

DPvT  -  ~- 
bT 

A=  Af 
E 


range  to  source  (km), 
vertical  correlation  length  (m), 

average  sound  speed  along  signal  path  of  propagation  (m/s), 

average  partial  derivative  of  sound  speed  with  respect  to  temperature  along  signal  path 
[(m/s)/°C], 

strength  of  temperature  field  (°C2/m),  and 
the  corresponding  environmental  factor. 


If  E  is  input  as  zero,  it  is  computed  from  values  of  VL,  CV,  DRT,  and  A.  If,  in  this  case,  any  of 
the  three  parameters  CV,  DRT,  or  A  are  input  as  zero,  those  equated  to  zero  are  determined  from  NP 
(input  variable)  profile  sets.  Each  profile  set  is  located  at  some  range,  RAN,  and  consists  of  one  or 
more  of  four  profiles:  sound  speed  (C),  temperature  (T),  Brunt-Vaisala  frequency  (BV),  and  salinity 
(SALT)  (°/oo)  versus  depth  (D).  For  each  of  the  NP  ranges  (RAN),  a  profile  set  is  input  by  reading 
NC  card-image  lines  of  the  form 

D(I),  C(I),  T(I),  BV(I),  SALT (I), 
where  NC  is  the  number  of  depth  points  and 


D  (I) 
C(I) 

TO) 

BV  (I) 
SALT(I) 


/th  depth  point  (m), 
sound  speed  at  D(I)  (m/s), 
temperature  at  D (I)  (°C), 

Brunt-Vaisala  frequency  x  103  at  D(I)  (10-3  rad/s), 
salinity  at  DO)  (°/oo) 


In  all  cases,  a  sound-speed  profile  (D(I)  and  C(I)  values]  must  be  specified  or  computed  for 
depth-averaging  purposes  (see  below).  However,  unless  CV  is  input  as  zero,  the  CO)  values  are  not 
used  to  determine  CV.  The  required  sound  speeds  C(I)  may  either  be  specified  or,  if  input  as  zero,  be 
computed  from  Wilson’s  equation  [32],  corresponding  input  temperature  TO)  and  salinity  SALTO) 
values,  and  a  computed  pressure  value  [33].  Zero  salinity  values  are  reset  to  35. 

If  CV  is  input  as  zero,  we  determine  the  average  sound  speed  by  first  averaging  the  CO)  values 
over  depth  at  each  range  then  averaging  the  individual  depth  averages  over  the  ranges  of  the  NP 
sound-speed  profiles. 

If  DRT  is  input  as  zero,  a  depth-range  average  temperature  T  is  similarly  found.  The  derivative 
DRT  -  bc/bT is  found  from  Wilson’s  equation,  i.e., 

DRT  -  4.623  -  0.1092T.  (31) 


If  A  is  input  as  zero,  a  depth-range  average  Brunt-Vaisfila  frequency  (xlO3)  BV  is  similarly  found 
and  the  field  strength  A  is  found  from  [4] 

A  -  BV  x  10-7.  (32) 

Similarly  to  the  treatment  of  sound-speed  profiles,  if  the  BV  (1)  values  of  any  profile  set  are  input  as 
zeroes  (and  A  is  input  as  zero),  they  ate  computed  from  corresponding  input  sound  speed  CO),  tem¬ 
perature  TO),  salinity  SALTO),  and  computed  density  values  [33], 
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The  depth  averaging,  performed  at  each  range,  is  an  approximation  to  averaging  over  the  cyclic 
path  of  a  dominant,  horizontal  signal  arrival.  A  receiver  depth,  RD,  is  specified  for  which  the  sound 
speed  is  determined  from  the  input  sound-speed  profile  [i.e.,  C(I)  values  -  recall  that  C(I)  must  be 
specified].  Next,  a  reciprocal  depth  is  found.  If  a  true  reciprocal  de,Hh  dees  not  exist,  the  reciprocal 
depth  is  taken  to  be  either  the  ocean  surface  or  the  ocean  bottom,  whichever  is  appropriate.  The 
reciprocal  depth  is  assumed  to  be  opposite  in  depth  from  the  minimum  sound  speed  (durt  axis).  That 
is,  if  the  receiver  depth  is  above  the  depth  of  the  minimum  sound  speed  (duct  axis),  the  reciprocal 
depth  is  assumed  to  be  below  the  duct  axis,  and  vice  versa. 

The  depth  averaging  described  above  assumes  a  horizontal  arrival.  To  account  for  a  dominant 
arrival  that  does  not  arrive  horizontally,  the  receiver  depth  RD  need  only  be  leinte'pretcd  as  a  depth 
where  the  dominant  ray  path  becomes  horizontal  (vertex  depth). 

The  foregoing  discussions  assumed  that  the  environmental  parameter  E  is  input  as  zero,  and  also 
that  at  least  one  of  the  parameters  CV,  DRT,  and  A  is  input  as  zero.  If  a  positive  value  is  input  for  E, 
then  VL,  DRT,  and  A  are  ignored.  However,  a  value  for  CV  is  required  to  convert  from  frequency  to 
wavenumber.  If  zero  is  input  for  CV,  its  value  is  determined  from  input  profile  sets  in  an  identical 
manner  as  in  the  case  when  E  is  input  as  zero. 

Subroutine  SPACOH 


This  subroutine  evaluates  the  NRL  coherence  model  using  the  environmental  parameters  deter¬ 
mined  by  subroutine  MEDIUM.  Subroutine  SPACOH  evaluates  the  function 

C(s)  -  exp  l—  Ek^xs^2]  (33) 

at  4000  equally  spaced  samples  of  the  separation  distance  s,  between  s  —  0  and  su,  where 

C(s«)  -  exp  (-32)  -  1.26  x  lO"14.  (3-  : 

The  value  exp(-32)  was  chosen  to  reduce  the  truncation  error  in  approximating  the  infinite  Fourier 
transform  of  C(s)  by  a  finite  transform.  The  large  number  of  spatial  samples  reduces  the  aliasing 
effects  of  approximating  a  continuous  transform  by  a  discrete  transform. 

A  subset  of  the  sampled  coherence  function  is  printed.  As  an  option,  a  CALCOMP  plot  of  the 
spatial-coherence  function  is  generated.  Also,  as  an  option,  the  sample  values  are  written  onto  an  out¬ 
put  file  on  unit  LCOH.  This  file  may  be  used  as  an  input  file  for  later  reruns  of  program  COHORT. 

Subroutine  WAVCOH 


This  subroutine  estimates  the  signal  (coherence)  spectrum  by  evaluating  the  discrete  Fourier 
transform  of  sampled  values  of  the  spatial-coherence  function.  These  sample  values  are  either  those 
found  in  subroutine  SPACOH  or  comparable  values  read  from  an  input  file. 


The  signal  spectrum  S(k )  is  evaluated  at  Ac  —  0  and  at  symmetrical,  equally-spaced  samples  span¬ 
ning  the  domain  -lc0  <  k  <  k0.  The  sample  spacing  in  k- space  is 


Aj  " 


JL 


(35) 


where  A  is  the  separation  in  k-spacc  of  independent  beams  of  an  uniformly  weighted  line  array  with  N 
elements  equally  spaced  a  distance  d  apart.  It  also  corresponds  to  the  distance  between  pattern  nulls  of 
such  an  array.  The  parameter  Nf  (called  a  fill  factor)  allows  for  a  finer  sampling  than  A.  Sample  spac¬ 
ing  is  discussed  further  in  the  next  subsection,  which  describes  the  array-performance  subroutines. 
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» < 


The  sample  spacing  A,  is  determined  from  the  input  parameters 

NPGS,NFILL,DISTPG,SS, 

where 

NPGS  =  N  =  number  of  array  elements,  or  group  modules; 

D1STPG  —  distance  (m)  between  array  elements  or  group  centers; 

SS  =  c  =  sound-speed  (m/s)  at  the  array;  and 
NF1LL  -  Nf  =  fill  factor. 

The  total  number  of  samples  must  not  exceed  4000. 

The  signal  spectrum  is  normalized  by  S(0)  and  referred  to  as  the  field  pattern.  Sample  values  of 
the  field  pattern  are  printed.  As  options,  the  field  pattern  can  be  plotted  and/or  written  to  an  output 
file  on  unit  KCOH. 


If  a  comparable  signal  spectrum  was  computed  by  a  previous  computer  run  and  written  to  an  out¬ 
put  file,  that  file  can  be  read  by  subroutine  WAVCOH,  as  an  option,  and  the  spectrum  printed  and, 
optionally,  plotted.  This  option  is  an  alternative  to  computing  the  spectrum  from  either  input  samples 
of  the  spatial  coherence  function  or  from  one  computed  from  input  environmental  parameters. 

Array-Performance  Subroutines 

Array  performance  can  be  predicted  for  any  of  several  array  configurations.  The  receiver  is 
assumed  to  be  an  array  of  identical  hydrophone  groups.  A  hydrophone  group  may  consist  of  a  single 
omnidirectional  hydrophone,  or  it  may  be  in  itself  a  linear  array  of  equally  spaced  hydrophones  which 
cannot  be  steered  from  broadside.  Nominally,  the  receiver  is  assumed  to  be  a  horizontal  line  array  of  N 
identical  hydrophone  groups  equally  spaced  a  distance  d  apart.  The  performance  of  other  array 
configurations  (e.g.,  distorted  line  array,  nonequal  spacing)  can  be  predicted  by  the  careful  specification 
of  iV  and  d  (to  determine  sample  spacings  and  number  of  beams,  see  below)  and  by  reading  in  the  array 
pattern.  In  all  cases,  the  receiver  pattern  is  the  product  pattern  of  the  array  pattern  times  the  group  pat¬ 
tern ,  as  shown  in  the  section  "Beam  Patterns." 


The  field  pattern  (signal  spectrum)  computed  in  subroutine  WAVCOH  serves  as  an  input  to  the 
array-performance  subroutines.  The  beam  patterns  and  the  field  pattern  must  be  specified  at  the  same 
sampled  k  values:  A  —  0  and  symmetrical,  equally  spaced  values  spanning  -k0  <  k  <  k0.  The  sample 
spacing  is  gi'  cn  by 


A,  - 


(36) 


which  is  determined  in  subroutine  WAVCOH  from  input  parameters.  The  distance  A  is  the  spacing  in 
*-space  of  the  independent  beams  of  a  uniformly  weighted  line  array  of  equally  spaced  elements.  The 
output  from  N  such  beams  is  produced  by  the  application  of  a  discrete  Fourier  transform  to  element 
data  (this  is  called  FFT  beamforming).  For  the  broadside  beam  pattern  of  such  an  array,  the  spacing 
would  produce  pattern  samples  only  at  the  main-beam  axis  and  at  nulls.  The  fill  factor  Nf  is  used  to 
produce  a  finer  sampling.  There  arc  two  subroutines  concerned  with  array  performance  prediction: 
BEAMS  and  CONVLV. 


Subroutine  BEAMS 

This  subroutine  computes,  prints,  and.  optionally,  plots  the  array  and  group  patterns.  The  input 
parameters  of  subroutine  BEAMS  are 
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NPAT,  NPFILL,  NPH,  and  DISTPH. 


The  array  pattern,  Ba(k),  may  either  be  read  from  catd-image  data  or  be  computed.  If  the  input 
parameter  NPAT  <  0,  the  parameter  NPFILL  is  ignored  and  the  array  pattern  is  computed  from 

Ba(k)  “  5aU°  Sin  }  "  N  sin  (kd/2) 

which  is  the  theoretical  pattern  for  a  perfectly  linear  array  of  equally  spaced  hydrophones. 


(37) 


If  the  input  parameter  NPAT  >  0,  the  parameter  NPFILL  is  required  to  be  identical  to  NFILL 
(Nf)  as  a  safety  check,  and  the  array  pattern  is  input  from  card-image  data  in  the  following  manner: 


•  NPAT  values  of  the  array  pattern  are  read. 


•  Pattern  values  are  read  in  order  of  increasing  angles  whose  sines  are  integer  multiples  of 


A, 

*o 


_1_ 

Nf 


X 

Nd 


•  The  NCENth  sample  corresponds  to  0  -  0,  where  NCEN  —  NPAT/2  if  NPAT  is  even 
and  NCEN  -  (NPAT  +  0/2  if  NPAT  is  odd. 


If  the  input  array  pattern  does  not  span  — tt/2  to  +ir/ 2,  it  will  be  repeated  (in  A-space)  in  order 
that  the  array  pattern  will  be  specified  at  the  same  sample  points  as  the  field  pattern.  Note  that  the  pat¬ 
tern  is  strictly  repeated.  For  example,  if  the  first  and  last  (NPAT)  pattern  values  are  zero  (nulls),  then 
two  successive  zero  values  will  be  generated  whenever  the  pattern  is  repeated. 


The  group  pattern,  BAk ),  is  computed  from 


B,(k)  - 


sin  (Mkde/2) 


M  sin  (kdJ2) 


where 

M  -  NPH  «*  number  of  hydrophones  per  group,  and 
</,  -  DISTPH  —  group  hydrophone  separation  (m). 


(38) 


Subroutine  CONVLV 


Subroutine  CONVLV  convolves  the  field  pattern  (signal  spectrum)  with  the  array-group  product 
pattern  of  the  receiver  to  calculate  the  loss  of  signal  gain  (re  ideal)  for  each  of  Abeam  directions.  The 
N  beam  directions,  or  angles,  bracket  the  signal  direction  and  arc  spaced  at  intervals  of 


S  - 


(39) 


in  k  space.  The  N  beams  arc  numbered  sequentially  from  negative  to  positive  such  that  the  main  axis 
of  beam  number  0  coincides  with  the  signal  direction.  The  convolution  can  be  performed  for  several 
signal  directions. 


The  input  parameters  of  subroutine  CONVLV  arc 


NSHF  (l).  I  -  1.  NCENT. 
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which  specify  NCENT  signal-arrival  directions,  for  each  of  which  the  field  and  product  patterns  are  con¬ 
volved.  For  each  signal  case,  the  signal-arrival  angle  corresponds  to  the  main-axis  angle  of  the  beam 
whose  beam  number  is  given  by  NSHF(l),  i.e., 

sine,  -  NSHF(I)  •  (40a) 


or  in  wavenumber  space. 


k. 


NSHF(I) 


kp\ 

~Nd‘ 


(40b) 


For  each  signal  case,  the  (signal)  field  pattern  is  also  shifted  from  broadside  to  the  arrival  direction  ks, 
as  shown  in  Fig.  7.  The  receiver  pattern  is  steered  sequentially  through  the  N  beam  angles  bracketing 
the  signal  direction.  For  each  steerage,  the  shifted  field  and  product  patterns  are  multiplied  and  scaled 
to  estimate  the  loss  in  signal  gain  for  that  beam.  This  process  is  repeated  for  subsequent  signal  arrival 
directions  if  NCENT  >  1. 


~kn  0  ka 


Fig.  7  —  Field  pattern  for  arrival  direction  k, 


Recall  that  only  the  array  pattern  and  not  the  group  pattern  can  be  steered,  and  that  the  product 
pattern  is  normalized  by  the  group-pattern  response  in  the  steered  direction  (to  maintain  a  constant 
main-axis  response).  For  example,  when  the  array  is  steered  to  the  direction  k„  the  receiver  product 
pattern  is 

B,(k.  k ,)  -  Bt(k  -  kf)Bt(k)/B,(k,),  (41a) 

or 

B,(k,.  k t)  -  B,( 0).  (41b) 

For  each  signal-arrival  direction  k,,  subroutine  CONVLV  prints  the  parameters 

S,  —  beam  angle  (re  signal  arrival-angle). 

Llk„  k,)  —  loss  of  signal  gain  (re  ideal), 

B,(k,)  —  product  pattern  normalization. 

fij(0)  “  main-axis  response  of  array  pattern,  and 
S(k,  -  k,)  —  signal-field  strength  at  beam  angle  9, 

as  functions  of  beam  number  l  Also,  as  an  option,  a  plot  of  loss  versus  (relative)  beam-angle  is  gen¬ 
erated. 

After  all  NCENT  signal  cases  are  processed,  the  tosses  at  each  f relative)  beam-angle  arc  averaged 
(over  she  signal  cases  or  array  shifts).  These  averages  are  printed  and.  cptionally.  plotted. 
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OPERATING  INSTRUCTIONS 

This  section  describes  the  card-image  input-data  structure  of  program  COHORT  and  presents  a 
sample  set  of  JSL  instructions  for  executing  the  program  on  the  ASC  resident  at  NRL. 

Card-Image  Input 

This  section  begins  with  a  card-by-card  description  of  the  COHORT  card-image  input  data  fol¬ 
lowed  by  a  summary  of  this  structure.  Next  is  a  table  of  the  input  structures  for  the  three  basic  options 
of  program  COHORT.  A  flowchart  outlining  the  various  program  options  and  a  sample  set  of  JSL 
instructions  conclude  the  section. 

Card-Image  Data 

The  following  describes  the  input-data  structure  for  running  program  COHORT: 

Card  1  inputs  are  KBEAM,  LCOH,  KCOH,  AND  KPLOT,  with  format  (415),  where 

KBEAM  -  code  for  selecting  program  options, 

LCOH  *  logical  unit  for  reading  or  writing  spatial  coherence, 

KCOH  -  logical  unit  for  reading  or  writing  coherence  spectrum,  and 

KPLOT  species  options  for  plotting  coherence. 

Options  are  as  follows: 


•  KBEAM  —  l,  compute  spatial  coherence  and  determine  beam  outputs, 

LCOH  >  0,  write  spatial  coherence  onto  unit  LCOH; 

KCOH  >  0,  write  coherence  spectrum  onto  unit  KCOH; 


•  KBEAM  -  2,  input  coherence  and  determine  beam  outputs. 


LCOH  >  0.  read  spatial  coherence  from  unit  LCOH  and  compute  coherence 
spectrum,  if  also  KCOH  >  0,  write  spectrum  onto  unit  KCOH, 
LCOH  5s  0.  read  coherence  spectrum  from  unit  KCOH  (must  tie  positive). 


•  KBEAM  ->  3.  compute  coherence  only. 


LCOH  >  0.  write  spatial  coherence  onto  unit  LCOH. 

KCOH  <  0.  compute  spectrum,  but  do  not  write  to  an  output  file, 
KCOH  —  0.  do  not  compute  spectrum. 

KCOH  >  0.  compute  spectrum  and  write  onto  unit  KCOH; 

•  KPLOT  ”  0.  do  not  generate  pk>u  of  coherence. 

•  KPLOT  —  1.  plot  only  spatial  coherence. 


•  KPLOT  ~  2.  plot  both  spatial  coherence  and  its  spectrum. 


•  KPLOT  —  3.  plot  only  spectrum 


Card  2  input  is  P.  with  formal  (E10  2).  where 
F  —  acoustic  frequency  (Hr) 
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Card  3  inputs  are  R,  VL,  CV,  DRT,  A,  and  E,  with  format  (4F10.2,  2E10.2). 

These  parameters  specify  the  NRL  coherence  model,  where 

R  =  range  (km), 

VL  =  vertical  correlation  length  (m), 

CV  =  average  sound  speed  (m/s), 

DRT  —  average  partial  derivative  of  sound  speed  with  respect  to  temperature 

[(m/s)/°C], 

A  =»  strength  of  temperature  fluctuations  (°C2/m),  and 

E  =  corresponding  environmental  parameter. 

NOTE:  This  card  is  read  only  if  KBEAM  =-  1  or  3  (coherence 

functions  are  read  in  rather  than  computed  when  KBEAM  **  2). 

Card  4a  inputs  are  NP  and  RD,  with  format  (15,  F10.2),  where 

NP  —  number  of  profile  sets  to  follow,  and 
RD  —  receiver  depth  (m). 

Card  4b  inputs  are  NC,  RAN,  and  FLAT,  with  format  (15,  2F10.2),  where 

NC  —  number  of  depth  points  to  follow, 

R  AN  -  range  of  profile  set  (km),  and 

FLAT  latitude  of  profile  set  (deg)— ignored  unless  sound  speeds  are  calculated. 

Card  4c  inputs  are  (D(I),  S(I),  T(l),  BV(I),  SALT(l),  I  -  1,  NC),  with  format  (5F10.2,/),  where 

D(I)  =  Ith  depth  (m), 

S(I)  ”  1th  sound-speed  (m/s), 

T(I)  *  Ith  temperature  (°C), 

BV(1;  =  Ith  Brunt-Vaisula  frequency  (r/s  x  103),  and 

SALT (I)  —  Ith  salinity  (°/oo)— defaults  to  35. 

NOTE:  One  Card  4b  followed  by  NC  Cards  4c  are  read  for  each  of  NP  ranges. 

Cards  4a  through  4c  are  read  only  if  KBEAM  —  1  or  3 
if  E  =  0  and  any  one  of  CV,  DRT,  and  A  is  zero  or  if 
E  >  0  and  CV  -  0. 

See  the  section  on  Subroutine  MEDIUM  for  a  complete 
description  of  options  for  specifying  parameters  of  the 
coherence  model. 

Card  5  inputs  are  NGPS,  NFILL,  DISTPG,  and  SS,  with  format  (215,  2F10.5),  where 

NGPS  —  number  of  group  modules, 

NFILL  —  fill  factor  for  evaluating  coherence  spectrum, 

DISTPG  =  distance  between  center  points  of  any  two  acfjacent  group  modules  (m),  and 
SS  —  sound  speed  at  the  array  (m/s). 
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NOTE: 


Sample  coherence-spectrum  values  are  equally  spaced  in  *  space  b, 

4, - CWF__  w  , 

NF1LL  *NGPS*DISTPG  NfNd  = 

This  card  is  read  only  if 

"  '  (spi‘l1,ial  coherence  and  spectrum  computed),  or 
compmed)'  or  W°H  "  °  (spaltal  «*  «ctrum 

®°  computed)"  3  a"d  KC°H  "  °  (spa,ial  “h«"“  ««> 


If  KBEAM  -  3,  no  more  cards  follow. 


FI0.2K  whe,ei"PU'S  "  NPAT'  »»,  NCENT,  and  DISTPH,  with  forma,  (515, 


npat 

NPAT 

NPFILL 


NPH 

NCENT 

DISTPH 


’  * - * 


0,  read  in  number  of  array-pattern  values; 

0,  compute  standard  array  pattern; 
fill  factor  of  sample  beam-pattern  values; 
ignore  if  NP AT  <  0;  if  NPAT  >  0  NPFILL 
must  equal  NFILL,  which  is  checked  for 
normally,  NPAT  -  NGPS’NPFILL; 
number  of  hydrophones  per  group  module; 
number  of  beam  shifts  desired  (0  <  NCENT  <  9)-  and 
distance  between  individual  group  phones  (m).  ’ 


Card  7  inputs  are  (NSHF(I),  I  -  ,,  NCENT),  with  forma,  (1615),  where 

NSHF(I)  -  amount  of  shift  of  Ilh  beam  shift  in  units  of  the  spacing  of  the 
independent  array-pattern  beams  in  it  space-  this 
spacing  is  given  by 

VSS/F  *oA 
NCPSWSTPd’  Nd " 


A  = 


Card  8  input  are  (PAT(I),  1  =  1,  NPAT),  with  format  (8E10.2),  where 
PAT(I)  =  array  pattern  in  real  units  (i.e.,  not  dB). 

NOTE:  Array  pattern  values  must  be  equally  spaced  in  k  space  by  A, 

This  card  is  read  only  if  NPAT  >  0. 

Card  9  input  is  KASPLT,  with  forma.  (15),  where 


KASPLT 

KASPLT 

KASPLT 

KASPLT 

KASPLT 


specifies  options  for  plotting  signal  loss; 

0,  no  plots  generated; 

L  plot  array  signal  loss  for  each  of  NCENT  signal  cases 

t’  p;ot  loss  for  each  casa  and  average  loss  for  all  cases;  and 
a,  plot  average  signal  loss  for  totality  of  losses. 
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The  following  is  a  summary  of  the  card-image  data  read  by  program  COHORT: 


Card 

Inputs 

Format 

1 

KBEAM,  LCOH.  KCOH,  KPLOT 

(415) 

2 

F 

(F10.2) 

3 

R,  VL,  CV,  DRT,  A,  E 

(4F10.2,  2E10.2) 

4a 

NP,  RD 

> 

(15,  F10.2) 

4b 

NC,  RAN,  FLAT 

■  repeat 

(15,  2F10.2) 

4c 

D(I),  S(I),  T(I),  BV(I),  SALT(I) 
(repeat  NC  times)  t 

NP  times 

(5F10.2,/) 

5 

NPGPS,  NFILL,  DISTPG,  SS 

(215,  2F10.5) 

6 

KBMPLT,  NPAT,  NPF1LL,  NPH,  NCENT,  DISTPH 

(515,  F10.2) 

7 

(NSHF(I),  1  =  1,  NCENT) 

(1615) 

8 

(PAT(I),  I  =  1,  NPAT) 

(repeat  NPAT/8  times) 

(8E10.2) 

9 

KASPLT 

(15) 

The  following  is  a  table  listing  the  required  and  optional  input-data  cards  for  the  program  options 
determined  by  the  input  value  of  KBEAM. 

Option  Table 


The  following  table  shows  the  input  structures  for  the  three  basic  options  of  program  COHORT: 


Card 

KBEAM  =  1 

KBEAM  -  2 

_ 

KBEAM  =  3 

— 

Subroutine 

Required 

Optional 

Required 

Optional 

Required 

Optional 

1 

X 

X 

X 

— 

2 

X 

X 

X 

— 

3 

X 

X 

MEDIUM 

4a 

X 

X 

MEDIUM 

4b 

X 

X 

MEDIUM 

4c 

X 

X 

MEDIUM 

5 

X 

X 

X 

WAVCOH 

6 

X 

X 

BEAMS 

7 

X 

X 

CONVLV 

8 

X 

X 

BEAMS 

9.__ 

X 

X 

- 

Flow  hart 

The  following  flowchart.  Fig.  8,  relates  the  card-image  data  with  the  program  options  prescribed 
by  the  parameters  KBEAM,  LCOH,  and  KCOH. 
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Fig.  8  -  Flowchart  for  COHORT 
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Fig.  8  (Continued)  —  Flowchart  for  COHORT 
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Fig.  8  (Continued)  -  Flowchart  for  COHORT 
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READ  CARD  #9:  KASPLT 


KASPLT  3  2  OR  3 
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Sample  JSL  Commands 


The  following  is  a  sample  set  of  JSL  instructions  for  executing  COHORT  on  the  ASC  system: 

/  JOB,  jobname,  account  number,  usercode,  CAT  =  1,  OPT  —  (T,  R) 

/  LIMIT  BAND  =  SO,  SEC  =  120 
/  JSLOPTS  OPT  =  (F,  L,  M) 

/  PD  YOU,  USERCAT/D—/B— /usercode/ node 

/  ASG  SYS.OMOD,  USERCAT/D81/L60/COHORT/OBJECT.  USE  =  SHR 
/  ASG  FT01F001,  YOU/SCOH,  USE  -  SHR* 

/  ASG  FT02F001,  YOU/KSPEC,  USE  -  SHR* 

/  FD  FT06F001,  BAND  -  2/20/2 
/  LNK 

/  FXQT  CPTIME  -  10000 


card-image  data 


/  FOSYS  FT06F001 

/  CAT  YOU/SCOH,  ACNM  -  FTOlFOOlt 
/  CAT  YOU/KSPEC,  ACNM  -  FT02F001t 
/  FOSYS  FT59F001,  TYPE-PLOT,  FORM-00* 
/  EOJ 


*t»The  inclusion  of  these  commends  depends  upon  user -selectable  options. 

♦If  it  is  specified  that  an  input  file  is  to  read  from  logical  unit  N.  then  a  JSL  command  of  the  form 
/  ASO  FTONFOOI,  YOU/PtLENAME.  USE-SHR 

must  be  supplied. 

♦  If  it  is  specified  that  an  output  file  is  to  be  written  onto  logical  unit  N,  then  a  JSL  command  of  the  form 
/  CAT  YOU/FILENAME.  ACNM  -  FTONFOOI 

must  be  specified 

Ml' one  or  more  CALCOMP  plots  arc  requested,  then  this  JSL  command  must  be  furnished. 
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Appendix 
SAMPLE  CASE 

This  appendix  describes  a  sample  execution  of  program  COHORT.  The  purpose  of  this  example 
is  to  provide  a  numerical  illustration  of  the  type  of  output  program  COHORT  is  capable  of  producing 
and  to  serve  as  a  guide  to  the  user  initially  attempting  to  execute  the  program.  The  option  selected 
(KBEAM  — 1)  computes  both  signal  coherence  and  array  performance,  thus  utilizing  all  the  subroutines 
(or  functions)  of  program  COHORT. 

INPUT  DATA 


The  card-image  input  data  used  to  generate  the  sample  results  are  listed  in  Fig.  Al. 
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Fig.  At  —  Cird-imige  input  diu  for  simple  case 
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ENVIRONMENTAL  PARAMETERS 


To  compute  spatial  coherence,  the  signc.  source  was  assumed  to  be  located  at  a  range  of  2000  km 
and  to  transmit  at  an  acoustic  frequency  of  300  Hz.  To  compute  the  environmental  parameter 


E  —  1.7 


I  ML 

c  BT 


it  was  assumed  that 


and 


c-  1501.1  (m/s), 


ML 

d  T 


-  4.26  (m/s)/°C, 


/v  -  30  m; 

but  the  strength  of  the  random  temperature  fluctuations  A$  was  assumed  to  be  unknown.  To  estimate 
A $  (which  is,  effectively,  the  average  Brunt- Vaisala  frequency  along  the  propagation  path)  a  single 
profile  set  was  input.  The  profile  set  specified  sound  speed,  temperature,  and  salinity  at  each  of  44 
depths.  Recall  that  zero  sound  speeds  are  recalculated  from  temperature,  depth,  salinity,  and  latitude. 
The  profile-set  range  was  arbitrarily  set  to  zero  and  the  latitude  taken  to  be  32.5°  (ignored  here). 


The  receiver  depth  was  set  at  zero  (the  sea  surface),  which  resulted  in  the  Brunt-Vaisala 
frequency’s  being  averaged  between  the  sea  surface  and  a  depth  of  4418  m  (the  reciprocal  depth).  The 
computed  strength  parameter  was 

A }  -  1.48  x  10'7, 

which  resulted  in 


SIGNAL  COHERENCE 


E  -  6.09  x  10'". 


Next,  the  normalized  coherence  function  was  evaluated  as  a  function  of  horizontal  separation. 
The  optional  plot  of  the  normalized  spatial  coherence  function  is  shown  in  Fig.  A2.  Finally,  the  option 
to  compute  and  plot  the  wavenumber  spectrum  of  signal  coherence  was  exercised,  and  the  result  is 
shown  as  Fig.  A3. 


Fi*  A2  —  Normatucsl  sjutwl  adherence 
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Fig.  A3  -  Normalized  wavenumber  spectrum  of  signal  coherence 


BEAM  PATTERNS 

The  horizontal  receiving  array  was  assumed  to  consist  of  128  hydrophone  groups  equallv  spaced 
5  m  apart  (center-to-center)  Each  hydrophone  group  was  assumed  to  consist  of  four  omnidirectional 
hydrophones  spread  1.25  m  apart.  Corresponding  plots  of  the  array  and  group  beam-patterns  arc  shown 
in  Figs.  A4  and  A5.  respectively. 


f  tp  A*  -  A ttai  hcam 
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For  ihe  above  array  configuration.  independent  beams  (and  also  nulls)  are  equally  spaced  in 
wavenumber  space  by 

A  =  16it  x  1(T4  m'1. 

If  the  array  were  not  more  finely  sampled,  u  woulo  be  specified  at  the  main-beam  axis  (unit  response) 
and  at  nulls  'zero  response).  To  affect  a  meaningful  convolution,  the  sampling  rate  was  increased  by  a 
factor  of  2  (i.c.,  NULL ~2)  Array,  group,  and  signal  patterns  wete  sampled  (internally)  at  this  spac¬ 
ing 

ARRAY  PERFORMANCE 

Array  performance  was  predicted  for  each  of  three  signal-arrival  angles  As  required,  each  arrival 
angle  corresponded  to  the  mam-beam  aits  of  an  independent  beam  angle  of  she  array  pattern  (with  the 
0'  broadside  beam  angle  .serving  as  a  reference )  Recall  thal  these  beam  angles  arc  on  the  main-beam 
axes  of  spatially  independent  beams  that  would  result  if  the  array  were  perfectly  linear  'wh.ch  in  this 
case,  it  ts)  and  the  hydrophone  groups  each  consisted  of  a  single  omnidirectional  hydrophone  For  the 
example  considered  of  a  l2S-e!emcnt  array,  there  arc  12S  such  beams,  numbered  here  between  -63 
and  +64  with  beam  0  representing  the  broadside  beam  paitern 

The  three  signai-arnval  angles  selected  correspond  to  independent  beams  0  l broadside).  5.  and 
10  f  or  each  arrival  angle,  signal  gain  degradation  was  computed  for  each  of  128  independent  beams 
(directions!,  where  no  s  the  reference  beam  number  0  iv  directed  at  the  signal-arrival  direction  i  here- 
fore,  for  each  cave  the  signal-gam  degradation  is  computed  for  each  of  the  128  beam  angles  that  would 
result  from  phasc-stccrmg  the  array  to  the  signal  direction  and  performing  FIT  beamforming 

The  rcsulung  plots  of  signal-gain  degradation  for  the  three  signal  directions  arc  shown  in  Figs  A6 
through  AS  For  reference,  computed  signal-gam  degradations  arc  given  in  Table  A  l 


NORMALIZED  POWER  (dB) 
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Tabie  A 1— Signal-Gain  Degradalion 


Receiver 

Beam 

Signal-Gain  Degradation  (dB) 

Si 

B 

an i 

V 

Signal 

Beam 

5 

Signal 

Beam 

10 

0 

2.49 

2.49 

8 

29.36 

29.27 

29.18 

16 

35.65 

35.48 

35.31 

Finally,  signal-gain  degradation  for  each  beam  is  averaged  over  the  signal-arrival  cases.  The  result 
is  plotted  in  Fig.  A9. 

AVERAGE  BEAM  OUTPUT 


-0.50  -0.25  0  00  0.25  0.50 

SINE  THETA 

Fig.  A9  —  Ave'age  normalized  array  signal  gain 
for  the  three  signal-arrival  cases 
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